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Recently an effective Lagrangian for the interactions of photons, Nambu-Goldstone bosons and 
superfiuid phonons in dense quark matter has been derived using anomaly matching arguments. In 
this paper we illuminate the nature of certain anomalous terms in this Lagrangian by an explicit 
microscopic calculation. We also generalize the corresponding construction to introduce the axion 
field. We derive an anomalous axion effective Lagrangian describing the interactions of axions with 
photons and superfiuid phonons in the dense matter background. This effective Lagrangian, among 
other things, implies that an axion current will be induced in the presence of magnetic field. We 
speculate that this current may be responsible for the explanation of neutron star kicks. 

PACS numbers: 



I. INTRODUCTION 

It has been recently realized^ that some very unusual effects may take place in dense matter QCD in the presence 
of topological defects and/or external magnetic field. It is known that at large baryon density many global symmetries 
of QCD are spontaneously broken0,Q. In particular, it is expected that the chiral symmetry will be spontaneously 
broken. This leads to appearance of pseudo-scalar Goldstone bosons, which are generated out of the vacuum by 
axial chiral currents. The spontaneous breaking of global symmetries also leads to existence of topological defects: 
domain walls and strings in dense QCD0,0|. In a recent paper it was shown that the effective Lagrangian for the 
interaction of Goldstone bosons with the electromagnetic field in the presence of chemical potential fi contains terms, 
which imply that topological defects such as axial vortices and domain walls in dense QCD carry electric current and 
magnetization respectively. The corresponding effective Lagrangian was derived in 1] in a formal way by treating 
the fcrmion chemical potential as the zeroth component of a fictitious vector gauge field and appealing to chiral 
anomalies induced by V^. 

The main goal of this paper is twofold. First, we observe that one of the terms in the effective Lagrangian derived 
in p] implies that flux tubes in dense quark matter carry an axial current. We investigate the microscopic origin of 
this current. Second, we generalize the derivation^] to include the axion field which, if exists, may play an important 
role in astrophysics and cosmology. 

This paper is organized as follows. In section II, we will be studying axial currents on flux tubes in dense matter. 
We will confirm the existence of these currents in three ways: a) by appealing to the effective Lagrangian derived in 
0, b) directly from the chiral anomaly due to the fictitious field V^, c) from a microscopic calculation. These three 
methods agree. 

In Section III, we derive the anomalous axion effective Lagrangian using previously established methods. The 
corresponding effective Lagrangian describes the interactions of axions with photons and superfiuid phonons in the 
dense matter background. We speculate regarding some phenomenological implications of the obtained results. 



'Electronic address: mmetlits@physics.ubc.ca 
^Electronic address: arz@physics.ubc.ca 



2 



II. AXIAL CURRENT ON MAGNETIC FLUX TUBES 
A. Short Overview 

The appearance of fermion zero modes on topological defects is a mathematically rich phenomenon with links to 
index theory of elliptic operators and quantum anomalies, see original papers^ and reviewQ. In particular, it is 
known0,0,Q that flux tubes in 2 + 1 dimensional QED possess fermion zero modes, whose existence can be proved 
using trace identities and 2-dimensional Euclidean chiral anomaly equations. In the context of 2 + 1 dimensional QED 
these zero modes lead to a degeneracy of the flux tube ground state and induction of a Chern-Simons term in the 
effective bosonic action. However, when we consider massless QED in 3 + 1 dimensions, the zero modes are free to 
move up or down the flux tube, depending on their chirality and polarization. We show in this paper, that for massless 
Dirac fermions in 3 + 1 dimensions at finite fermion chemical potential /j, and temperature T in the background of a 
magnetic flux tube, the zero modes generate an axial current of J = along the flux tube, where e is the fermion 

charge and $ is the total magnetic flux. The current is topological in nature as it depends only on the total flux and 
not on the particular details of distribution of the magnetic field. This result is exact (at least if the magnetic flux 
does not fluctuate) as the contribution of all, but the zero modes to the axial current vanishes. 

As already noted, the appearance of axial current on flux tubes can be derived by using a somewhat different 
aspect of chiral anomalies. Here the anomaly resides directly in 3 + 1 dimensions and appears when one thinks of 
the fermion chemical potential as the zeroth component of a fictitious vector gauge field V^. Such a gauge field, as 
well as the ordinary electromagnetic field A^, contributes to the anomalous non-conservation of axial current. If the 
chiral symmetry is spontaneously broken in the system under consideration (say in dense QCD) and a corresponding 
Goldstone boson r\ appears, one requires the effective Lagrangian for r\ to reproduce the axial current non-conservation. 
As will be shown below, the effective Lagrangian for 77, originally derived in pj, then implies the appearance of axial 
current on flux tubes, which agrees with the microscopic zero mode result. We also show that, alternatively, one can 
obtain this result starting directly from the axial current non-conservation equation, without appealing to the effective 
Lagrangian for the 77 boson. This later method, as well as the microscopic derivation, imply that axial current on flux 
tubes appears even if the chiral symmetry is not broken spontaneously. 

We note that although the presence of zero modes on flux tubes in QED and the symmetry of higher energy 
states, which makes the microscopic calculations of this paper possible, have been previously knownQ, the problem 
of calculation of axial current on flux tubes at finite chemical potential and temperature, to our knowledge, has not 
been considered before. Moreover, the computation of the axial current at finite chemical potential using fictitious 
anomalies in 3 + 1 dimensions is certainly a new trick. This work confirms the validity of this trick byan explicit 
microscopic calculation, thus, supporting the original derivation of the anomalous effective Lagrangian|l|. We note 
the similarity of computation of axial current on flux tubes in QED undertaken in the present paper to the recent 
computation of electric current on cosmic strings at finite chemical potential^, which was also motivated by the 
anomalous terms in the effective Lagrangian for Goldstones in dense QCD. Both calculations rely on the idea of zero 
modes with fixed quantum numbers running along a 2-dimensional topological object uniform in the 3rd direction 
and yield similar results for the current. 

B. Anomalous Terms 

Here we briefly review the anomaly based arguments of pj, which indicate that when an axial- like symmetry is 
spontaneously broken at finite fermion chemical potential, the effective Lagrangian for the corresponding Goldstone 
mode receives a correction, topological in nature. We show that in the presence of a background magnetic field, this 
correction leads to the appearance of axial current on magnetic flux tubes. 

Consider QCD at large baryon chemical potential. As is well known, such a system spontaneously breaks various 
global symmetries of QC.D0, y], leading to the existence of Goldstone bosons, whose number and form depends on 
the number of light (massless) quark flavours Nf. Here, we consider a general neutral Goldstone boson 77, 1 with the 
following transformation properties under one of the diagonal axial symmetries of QCD 2 : 

Va^e iQ »^Va, + (1) 



1 This need not be the "rf boson of QCD. 

2 As was shown in[j|, all instanton effects are suppressed in high-density QCD and at very large baryon chemical potential the explicit 
breaking of the formally anomalous symmetry U(1)a becomes very weak. 
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Here Q a denotes the flavor content of the Goldstone boson, a — l.JV/, and 77 is created out of the vacuum by the 
current, 

r = J2Q-^^ 5 ^ ( 2 ) 

a 

As is well known, it is useful to represent quark chemical potentials as the zeroth components of a fictitious vector 
gauge field — (1,0). Then the coupling of quarks to and to the usual electromagnetic gauge field takes the 
form: 

C = J2(^-e a A^ al ^ a (3) 

a 

where fi a and e a are quark chemical potentials and electromagnetic charges respectively. The anomaly equation for 
the current in the background of fields and takes the form: 

V - e^(C nAA F^F Xa + C nAV V^F Xa + C v vv V^V X<7 ) (4) 
where the field tensors Fx a , V\ a are defined as, F\ a = d\A a — d a A\, V\ a = d\ V a — d a V\, and the coefficients, 

C nAA = -N c J2j^r, C vAV = N c J2 6a ^? a > C nvv = -N c J2j^r (5) 

a a a 

The anomalous current non-conservation must be reproduced in the effective Lagrangian for the neutral Goldstone 
boson 77. Thus, as was shown inpj, the effective Lagrangian for 77 pics up the following anomalous term describing its 
interaction with the fields A^ and V^: 

L n = L° + 2d^ Xa {C 1lAA A y F x « + C vAV V„F Xa + C vVV V v V Xa ) (6) 

Here L° is the standard, non-anomalous part of the effective Lagrangian for 77, which transforms as L° — ► L° — j^d^Q 

under QJ. We now restore the fictitious field to its true value — (1,0). Then the last term in (JSJ) vanishes 3 , 
and we are left with, 

L n = L° v - C vAA r, eT^F^Fx* + 4C nAV Vr] ■ B (7) 

where B is the magnetic field. The first anomalous term in eq. JJJ describes the usual anomalous decay of a Goldstone 
boson to two photons, and is absent on the classical level, if there is no background electric field present. We now 
concentrate on the second anomalous term in JJJ, which does not occur in vacuum (at \x = 0): 

L n = L° + 4C vAV Vri ■ B (8) 

One effect of the new anomaly term originally discussed in is the magnetization of the domain walls formed by 
the 77 field (such domain walls are possible, say, if 77 is associated with spontaneous breaking of the U(1) A symmetry, 
which is also explicitly slightly broken by instantonsQ). Here we discuss a different consequence of this term. Let's 
vary the action obtained from Lagrangian JSJ), with respect to 7/ — » 77 + 9, to derive the classical equations of motion. 
By construction, L° — > L° — j^d^9, hence, 



5S = j d*x{-fd^e + AC riAV B % d l 9) = J (fx drf" 6 + j dt J dS t (-f + AC nAV B l )9 (9) 

Here the surface integral is over the boundary of the region R where our dense matter is realized. So, as V ■ B = 0, 
the anomalous term does not contribute to the equation of motion d^j^ = 0. However, if we do not restrict 9 to 
vanish on the boundary, we also obtain a boundary condition, 

J-dS = 4C nAV B-dS (10) 



However, as shown in ij this term can become important if the quark matter is rotating and/or superfiuid vortices appear. 
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Now, in the steady state situation, there is no build up of axial charge, and we have V • j = 0. Hence for any 
cross-section S of the region R let St be the part of dR such that dS = dSt,- Then, 

[ dS-j= j dS-j^4C vAV [ dS-B = AC nAV ! dS ■ B = iV c ^ ^"f" 3> (11) 

where $ is the total magnetic flux through the cross-section S. So we see that the anomalous term in eq. 10 implies 
the existence of an axial current flowing through the dense matter which is proportional to the magnetic flux. 

At this point we make the following important remark regarding the formula i|ll|) : the final expression for the 
current does not depend on the specific properties of the pseudo-Goldstone boson 77, such as its coupling constant f^. 
This is not due to our choice of units, and this is not a typo, so it leads us to weaken our assumption of spontaneous 
chiral symmetry breaking and existence of the 77 Goldstone. 

Indeed, the result 1)1 l[l can be derived in the following way for a generic system of massless fermions at finite 
chemical potential in the magnetic field, without appealing to the effective Lagrangian for the 77 Goldstone boson. 
Let's return to the anomaly equation |0J. We can think of the fictitious field Vq, as taking the value "1" inside the 
region R where the quark matter is realized and "0" outside (if we have a quark star, this is actually rather close 
to reality, since the interface between quark matter and vacuum is very narrow). Then, again, if no electric field is 
present, and if the axial charge density is time independent, eq. I@J takes the form, 

V-j = AC vAV W-{V B) (12) 

The right hand side of eq. I|12|) vanishes both inside and outside R, yielding V • j = 0. However, integrating 112(1 over 
a small Gaussian pillbox on the boundary of R, and recalling that in the vacuum outside R, j — 0, we obtain on the 
inner boundary of R, 

j-dS = AC vAV B-dS (13) 

This is the same result (|1U|) that we obtained by minimizing the action for the Goldstone 77. From this, we readily 
obtain the expression l|llfl for the total current across any cross-section of R. 

The last derivation does not use anywhere the spontaneous breaking of chiral symmetry and the existence of the 77 
Goldstone, and relies solely on chiral anomalies. Thus, when the dynamics of our problem are such that a Goldstone 
boson 77 exists, the appearance of axial current on flux tubes can be extracted from the effective Lagrangian © for the 
77 field. However, the existence of such axial current does not depend on spontaneous chiral symmetry breaking, but 
rather on the mere existence of chiral symmetry. This observation is confirmed microscopically in the next subsection, 
where the result (lll|) is reproduced in a simple QED-\ike system. 



C. Microscopic Arguments 

We will show in this section that the appearance of current on magnetic flux tubes at finite chemical potential derived 
in the previous section using the anomalous effective Lagrangian for 77, can be understood very simply microscopically 
within the following model. As we already saw at the end of the previous section, the existence of axial currents on 
flux tubes does not rely on spontaneous chiral symmetry breaking. Thus, we choose to ignore all the strong interaction 
effects leading to formation of the Goldstone field 77, and consider only the following Lagrangian, 

£ = iji{d^ + ieA^)Y^ - miptp + fMpj°ip. (14) 

which describes the interactions of a single light quark tp of mass m with a background electromagnetic field A^, at 
finite baryon chemical potential fi. Hence, the discussion in this section actually applies to any QED-\ike system at 
finite chemical potential. 

We are interested in the case of magnetic flux tubes, i.e. is static and the magnetic field B — V x A — B(x, y)z 
is uniform in the third direction z. Our goal is to compute the total axial current Jf = J d 2 x(ipj 3 ^/ 5 ip) along the flux 
tube. The Dirac Hamiltonian is, 

H = -i{di + ieA i ) 1 ° 1 l + mj° (15) 

and the Dirac equation becomes, 



-H R ip L + rmp R = Eijj L 
rmp L + H R ^ R = Eip R 



(16) 
(17) 
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where we use the conventions of Peskin and Schroeder and, 

H R = {-idi + eAi)^ 

So, ip L = ^(E - H R )ip R , where 

{H 2 R + m 2 )Vfl = E 2 ^ R 



(18) 



(19) 



Hence every eigenstate tpR of H R with eigenvalue e generates two solutions of the Dirac equation with energies 
E = ±Ve 2 + m 2 and, 



IpR 



2^-1 / ±((m 2 + e 2 )2 te)^ 



R 



((m 2 + e )2 ± e)zipi 



(20) 



Now we concentrate on the right sector HriJjr — etj)R. Due to invariance with respect to translation in z direction, 
we go to momentum eigenstates —id^R — P3,4>r. (we take the third direction to be periodic of length L, and take the 
limit L — > oo at the end of the calculation). In each momentum sector, the operator Hr takes form, 



Hr 
H 1 - 



p 3 a 3 + H 1 - 

(~id a + eA a )a a , a 



1,2. 



(21) 
(22) 



We note that {a 3 , H^} = 0. Hence, if |A) is a properly normalized eigenstate of H 1 - with eigenvalue A then ct 3 |A) 
is a properly normalized eigenstate of H with eigenvalue —A. So, all eigenstates of H 1 - with non-zero eigenvalues 
are of form |A), | — A) = er 3 |A), where A > 0. Also, u 3 maps zero eigenstates of H ± to zero eigenstates of H ± and 
hence we can classify all zero modes of H by their eigenvalue under c 3 . 

The eigenstates of Hr can now be expressed in terms of eigenstates of H . Clearly, [Hr,H^ ] = 0, so Hr only 
mixes states |A), | — A). For A > 0, we write, 



tpR = ci|A) +c 2 ct 3 |A) 



where c\,C2 satisfy: 



Hence e = ±-\/A 2 + p| and, 



A P3 
P3 -A 



(4(A 2 +^))" 



f ±sgn(p 3 )((X 2 +pl^±X)i 

\ ((a 2 +p!)^tA)^ 



(23) 



(24) 



(25) 



So each eigenstate of H 1 ^ with an eigenvalue A > generates two eigenstates of Hr. 4 
The zero modes of H 1 - are simultaneously eigenstates of Hr with eigenvalue, 

e = p 3 <r 3 



(26) 



Hence, when the mass m — > 0, zero modes of H 1 - become gapless modes of H capable of travelling up or down the flux 
tube depending on the sign of a 3 and on the chirality. We will shortly see, that at finite chemical potential, precisely 
these modes carry an axial current along the flux tube. 

The following quantity will be of particular importance to us: N — N + — 7V_, where N + and are the numbers of 
zero modes of iJ x with a 3 = 1 and a 3 = — 1 respectively. Observe, that if |A) is a zero mode of H 1 - with |A) = (it, v) 
then, 



Vv = 0, V^u = 



(27) 



4 There are known examplesQ, such as fermion number appearing on domain walls, when this is not strictly speaking true. Indeed, some 
of the energy levels of Hr are continuous rather than discreet and the correspondence discussed above between the eigenstates of H 1 - 
and Hr need not preserve the density of states. However, for the particular Hamiltonian Hr, it can be shown that if B(x) — > as 
x — * co sufficiently fast, this problem does not arise. 



6 



where, 

V = -id x - + e{Ax - %A 2 ) (28) 
Hence N + — dim(ker(T> : ' i )), N- = dim(ker(T>)), and, 

N = Index{H^) = N + -N_= dim(ker(V r )) - dim(ker(V)) (29) 



The index of the elliptic operator H 1 - has been computed in numerous works0, 0, |8j using two types of methods: 
i) complex analysis methods, ii) trace identities and axial Euclidean anomaly in 2 dimensions (this is particularly 
interesting in the light of our using 4 dimensional anomalies above to derive axial currents on flux tube at finite /i). 
The zero modes have also been computed exactly for some simple configurations of the gauge fieldQ. In general the 
index is given by: 

Index{H^) = %- (30) 

Z7T 



(31) 



Hence the index measures the number of flux quanta through the xy plane, which is in essence a topological quantity. 

Now let's proceed to compute the axial fermion current induced at finite chemical potential /i. For further generality, 
we also include the effects of non-zero temperature T . The axial current density in the third direction is given by, 

jl(x) = ?(i) 7 3 7V(i) = + ^ R ^ R (x) (32) 

We wish to compute the expectation value of the total current along the flux tube, J| = J d 2 x(j 3 (x)). At finite 
chemical potential and temperature we have, 

(j!(x)) = ^WWfeW = E(«((^ 2 + ™¥) + «(-(e 2 + (33) 

E e 

Here, n(E) — e ${E- 9 J)s^}(E) +1 is the usual Fermi-Dirac distribution, tpE are eigenstates of H with energy E, i^r^ are 
eigenstates of Hr with eigenvalue e, and we've used eq. (|20|l . The explicit form of ijjRe m terms of eigenstates of H 1 - 
implies, 



( J 5> = ^^^K(A 2 +^+m 2 )^)+n(-(A 2 +pl + m 2 )^)(^(A,p 3 )k 3 |^(A,p 3 )>+ (34) 

P3 A>0s=± 

+ ^EEW^3+^ 2 )^)+«(-(f3+™ 2 )^))<Ak 3 |A) (35) 

P3 A=0 

Here A > label eigenstates of H^, which generate eigenstates ^(A,^) of Hr with momentum p$ and eigen- 
value e± = ±\/A 2 + p 2 , while A = label the zero modes of H^. Now, let's evaluate the matrix element 
(ipji(X 7 p3)\a 3 \tp R (X, p^)} for A > 0. Using eq. l(2*3*j) and dropping the subscripts A, P3, s, we obtain, (i/)R\a 3 \ipR) — 
(|d | 2 + |c 2 | 2 )(A|cr 3 |A) + {c\c 2 + c x <%). Noting, (A|ct 3 |A) = (A| - A) = for A > 0, and using the explicit formula 
for ci, C2, we obtain, {'4'r\<j 3 \iPr) = sp3(X 2 + P3 2 )~^- This matrix element is odd in both P3 and s, hence the sum 
over all A > in eq. (|34|l vanishes , and only the zero modes of H 1 - contribute to J 3 . The zero modes carry a definite 
value of 173, so that (A|cr 3 |A) = a 3 . Thus, we are left with, 

Jf = (iY + -^)iEK(P3 + ™ 2 )^) + ^(-(P3 + ™ 2 )^)) = ^ m (T,M) (36) 

P3 

n m (T,») = J^(n((pl + m 2 )-i)+n(-(pi + m 2 )?)) (37) 

Here, n m (T,fi) is just the number density of one-dimensional two-component (Dirac) fermions of mass m at finite 
temperature T and chemical potential /1. Hence, our final result H36(l is topological in nature, since for each value of T 
and fx, it is sensitive only to the total magnetic flux and is independent of the particular distribution of the magnetic 
field. 
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Several limits of the result (|36(l are noteworthy. First of all, in the massless limit m — > 0, one has n(T, fi) — £ and, 

(38) 

If there are several species of quarks present, we can sum eq. (I38[l over quark flavours and colours to obtain the 
current J 3 of eq. J5J, which in the true dense QCD creates the 77 boson, 

J 3 = iV c V^#^$ (39) 

a 

This agrees with our result of the previous subsection, where we explicitly used the fact m = (and, hence, 
chiral symmetry) in assuming that the axial current conservation is violated only by anomalies. So, we see that the 
appearance of axial current on flux tubes, which was derived somewhat mysteriously in the previous section using the 
trick of fictitious chiral anomalies, is microscopically due to fermion zero modes. Our microscopic approach supports 
the validity of the fictitious chiral anomaly trick and serves as a check of the anomalous effective Lagrangian derived 
inQ. 

Let us note that the result Ij38|l is also independent of temperature for m = 0, which is a quite natural feature of 
a truly topological phenomenon. More explicitly, this fact is due to the special property of massless one-dimensional 
fermions, namely, their density at finite chemical potential is temperature independent. 

For arbitrary mass m ^ 0, the density of one-dimension fermions n(T, fi) is generally temperature dependent, so 
for simplicity we consider the limit T = 0. Then, n(0, fi) = \/ H 2 — m 2 /tt and, 



Ji= e -^^* (40) 

It is instructive to take the non-relativistic limit of eq. (|4UH . Writing, \x = m + fJ> nr , where the non-relativistic chemical 
potential /x nr <C m, 



jf « £Vp^$ (41) 

In the non-relativistic setting, Jf is just the spin S 3 , and for the case of uniform magnetic field, our result stems from 
the familiar fact that all Landau levels are doubly degenerate with respect to spin, except the lowest Landau level. It 
is amazing that this simple fact has such deep connections to chiral anomalies in 2 and 4 dimensions. 



III. AXION 



This section is devoted to the derivation of the anomalous effective lagrangian including the axion field using some 
previously developed methods. While the axion is considered to be one of the best dark matter candidates (see 
original papers[To | -|l3 | . and reviews |14|'). it has not been discovered yet. We derive novel low energy terms, which 
describe the interaction of the axion field with other light particles: photons and superfluid phonons in dense matter 
background. These terms may lead to phenomenologically important effects related to the axion astrophysics, which 
were not discussed previously. 

We define the 9 term in the fundamental QCD lagrangian in the standard way, Lg = -§2^G tiUa G ^„ . The existence 
of the term implies a violation of P, CP and T symmetries. However, there is no experimental evidence for P or CP 
violation in strong interactions. For example, CP violation in QCD would induce electric dipole moments of strongly 
interacting particles and there are stringent experimental limits on those quantities. Thus the absence of CP violating 
effects in QCD indicates a very small value for the parameter 0: why is so small? 

The most elegant resolution was proposed by Peccei and Quinn who assumed that the strong interactions Lagrangian 
has a global U(1)pq chiral symmetry |ldj|. Weinberg and Wilczek^lJ analysed the consequences of the Peccei-Quinn 

symmetry and noticed that the spontaneous breaking of a global chiral symmetry U(1)pq leads to a light pseu- 

2 

doscalar pseudo-Goldstone boson, called an axion, that will interact with topological charge density, -^^iG ilva G ^ va . 
In papers |l2lll3| two different types of the invisible axion models were suggested where it was demonstrated that the 
strong CP problem in QCD can be successfully solved with arbitrarily weak axion coupling constant. 

The only information which is relevant for us in what follows is the transformation properties of quarks under 
£7(1) pq chiral symmetry, 



(42) 
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where is the PQ charge for quark species ip a - We should note that leptons and Higgs bosons have also nontrivial 
transformation properties under U(1)pq symmetry, however this part is not essential for the present paper. We note 

2 

that the axion interacts with gluons through a triangle diagram, L agg = -j^^^G^ ua G llua . It also interacts with 

photons L a77 ~ ga-y^Q'F' 11 ' F^u with a coupling constant g a -yy ~ l//o Y] a Qa® e a expressed in terms of electric e a and 
PQ charges Q p ® of all quarks and leptons. In what follows, we identify the axion field with the dimensionless phase 
a{x) Physical, dimensional field a ~ f a a. 



A. Anomalous Axion Lagrangian 



The key observation here is as follows. The relevant for this work transformation properties of quarks under chiral 
rotations (JIJ and under PQ rotations (|42|l are very similar, Q a — ► . Therefore, we can literally follow our previous 
calculations (Section II) in order to derive the anomalous effective lagrangian for the axion field (which replaces the 
Goldstone field 77) in the presence of chemical potential /x. We use the same trickQ in this derivation by representing 
quark chemical potential as the zeroth component of a fictitious field V^. The result of this calculation is almost 
identical to ©, 

L a = L° a + 2d^ae^(C ajJ A u F x<J + C^ v V v F Xa + C aVV V„V x<J ) (43) 

where L° a describes all non-anomalous terms including the axion kinetic term /^(d^a) 2 , as well as different interaction 
terms of the axion with quarks and leptons. Coefficients Cq, 77 , C Q7 y, C a \?v can be easily extracted from the 
calculation of the triangle diagram and are given by, 

r \- zlQ P a Q r \- ^ a Q P ° r \- tiQ P a Q (M) 

6 «77 = -2^1^2-' Ca 1 V ~ 1^—^-2 ' W = "L l 67r 2 ' 

a a a 

where label a runs over all species (colours) of particles with nonzero PQ charges including quarks and leptons. The 
first term ~ C ai ^aF^ v in eq. (|43|l is the well known interaction between the axion and photons. 5 In particular, 
this term describes the axion decay into two photons. It also describes the axion<-> 7 transitions in the presence of 
the magnetic field B. The corresponding effect plays a crucial role in most axion search experiments. We shall not 
discuss this term in the present paper. 

Two other terms are new, and as far as we know, these terms have never been discussed in the literature. The 
last term in eq. (|43l) vanishes in the topologically trivial background. However, in the presence of superfluid vortices 
similar to the case discussed in pj , this term, among other things, describes the axion^superfluid phonon transitions. 
It might be phenomenologically important in rotating neutron stars. 



B. Possible Applications 

Let us concentrate on the middle term in brackets in eq. I|43|l . We rewrite it in terms of the physical fields in the 
following way, 

L a = 4C Q7 yV a . B (45) 

This interaction explicitly shows that a Peccei-Quinn current corresponding to the Peccei-Quinn symmetry will be 
induced in the presence of an external magnetic field. Indeed, one can literally follow the derivation to get the 
following expression for the Peccei-Quinn current, 

J?Q= f dS-j=4C aiV [ dS-B = J2 ea ^ Q ; Q <S> (46) 
Js Js o 27T 

where $ is the total magnetic flux through the cross-section S, and we have neglected the lepton chemical potentials. 
So we see that the anomalous term in eq. (|45|l implies the existence of a Peccei-Quinn current (|46|l flowing through 
the dense matter, which is proportional to the magnetic flux. 



5 In the literature devoted to the axion [T^l one typically includes in C a -n * ne effects due to the mixing of the axion with Goldstone fields 
such as the pion. We ignore this mixing for all qualitative discussions in what follows. 
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A few remarks are in order. First, as we noticed previously, the result l|46|l does not depend on f a similar to 
the previous case (|lip> when the result did not depend on f v . It may look very suspicious because one can make 
f a arbitrarily large, which corresponds to an arbitrarily small interaction of the original fermions with the axion. 
However, the point is that the current l|4()|) corresponds to the equilibrium state in an infinitely large bulk of matter. 
The relevant question for the present case is: what is the formation time for such a current? It is obvious that with 
f a increasing, the formation time also increases, such that there is no contradiction with eq. H46J) being independent 

Of fa- 

Our next remark is that the Peccei-Quinn i|4(j|) as well as the axial current are unique due to their topological 
nature. Indeed, even in the strongly interacting theory the axial current is persistent and non-dissipating. It 
means that even in such an unfriendly environment as the dense quark/nuclear matter in neutron stars the current 
does not dissipate due to re-scattering and can be effectively used to deliver information across the bulk of the star. 
Therefore, there is a unique opportunity here to use our topological currents (|46|) for delivering the asymmetry 
produced in the bulk of the star to solve the problem of neutron star kicks |15tll6| . 

The problem can be explained as follows. As is known, pulsars exhibit rapid proper motions characterized by a 
mean birth velocity of 450 ± 90 km/s. Their velocities range from 100 to 1600 km/s0, while their distribution 
leans toward the high- velocity end, with about 15% of all pulsars having speeds over 1000 km/s 16J. Pulsars are 
born in supernova explosions; therefore, it would be natural to look for an explanation in the internal dynamics of 
the supernova. However, three-dimensional numerical simulations |17j show that even the most extreme asymmetric 
explosions do not produce pulsar velocities greater than 200 km/s. Therefore, a different explanation should be 
found. The origin of these motions has been the subject of intense study and several possible explanations have 
been proposed. Many of the suggested mechanisms are capable ("in principle") to produce the required asymmetry. 
Indeed, in the presence of an external magnetic field, the produced neutrinos are automatically asymmetric with 
respect to the direction of B. However, the main common problem suffered by most suggested mechanisms is the 
difficulty of delivering the produced asymmetry to the surface of the star. Only in this case, the asymmetry may 
result in producing the proper motion of the entire star. To overcome the difficulty with delivery of the produced 
asymmetry to the surface, some proposals, for example, are based on new particles (such as a sterile neutrino) , which 
could escape from the bulk of the neutron star and deliver the asymmetry to the surface, see e.g. UM- 6 Our 
main observation here is that due to their topological nature, the currents l|ll|) . H46fl may be capable of delivering the 
required asymmetry (produced in the interior of the star) to the surface without dissipation|20|. 

IV. CONCLUSION. FUTURE DIRECTIONS. 

In this paper we have discussed the appearance of axial current on magnetic flux tubes at finite fermion chemical 
potential using several approaches. All of these approaches are weaved together by chiral anomalies. Microscopically, 
the current can be understood in terms of fermion zero modes on the flux tube. These fermion zero modes are in 
a certain sense themselves due to anomaly in 2 dimensional Euclidean field theory and have implications for the 
2 + 1 dimensional QED. Thus, we see that our trick with fictitious anomalies at finite chemical potential in 3 + 1 
dimensions, in a sense continues the propagation of anomaly from 2 to 3 to 4 dimensions. This is a common pattern 
in the study of anomalies. 

We anticipate a number of different applications of the derived anomalous effective low energy lagrangian for the 
Goldstone bosons and axion in dense matter. Some of them were mentioned in the original paperQ, others were 
mentioned in the present text. One specific application, which we believe deserves further study is the explanation 
of neutron star kicks [2 0|. In addition, novel anomalous effective lagrangian including the axion field might be quite 
important for analysis of a number of astrophysical problems, which would be the subject of a future work. 
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It is interesting to note that the interaction 1451 is a precise realization of the idea |lfll that the magnetic field may be correlated wit h 
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Dynamics, more precisely, anomaly, does the job of correlating the magnetic field with the weakly interacting axion current. 
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